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Dependence  of  Ocean  Heating  on  the  Distribution 
of  Spectral  Irradiance  in  the  North  Atlantic 

Robert  A.  Amone 
Gregory  E.  Terrie 
Paul  J.  Martin 

Naval  Research  Laboratory 
Stennis  Space  Center,  MS  39527 

ABSTRACT 

Subsurface  heatmg  rates  from  visible  solar  irradiance  were  computed  for  the  North  Atlantic  for  May  1979  during 
the  period  of  the  spring  bloom.  The  visible  spectrum  between  400  and  700  nm  accounts  for  a  substantial  fraction, 
about  43%,  of  the  total  solar  irradiance  at  the  sea  surface,  and  comprises  most  of  the  solar  irradiance  that 
penetrates  more  than  a  meter  into  the  sea.  The  mean  monthly  spectral  diffuse  attenuation  coefficients  and 
surface  solar  irradiance  were  used  to  compute  heating  rates  at  390,  440,  490,  540,  590,  640,  and  690  nm  over 
depth  inaements  of  0-5, 5-10, 10-15, 15-20,  and  20-35  m.  At  low  latitudes,  in  the  North  Atlantic,  significant  solar 
heating  occurs  at  depth  as  a  consequence  of  high  solar  irradiance  and  clear  waters.  In  the  northern  latitudes 
the  heating  is  confmed  near  the  surface  at  all  wavelengths  as  a  result  of  high  turbidity.  Significant  spatial  variation 
in  the  spectral  heating  rates  is  observed  as  a  result  of  chlorophyll  and  aerosol  patchiness. 


1.  INTRODUCTION 


Direct  solar  heating  of  the  upper  ocean  significantly  affects  physical  and  biological  processes.  The  distribution 
of  heat  near  the  oceans  surface  is  primarily  responsible  for  stabilizmg  the  surface  waters,  resulting  in  thermal 
and  biological  gradients  which  are  important  for  understanding  the  global  carbon  cycle  and  its  influence  on  global 
warming  K  The  rate  of  solar  heating  of  the  sea  is  controlled  by  the  incident  irradiance  and  the  optical  properties 
of  the  water.  The  variability  in  these  rates  is  responsible  for  the  ditumal,  seasonal,  and  climatological  scales  of 
thermal  evolution  in  the  upper  ocean  The  surface  heat  process  is  responsible  in  the  development  of  the  mixed 
layer.  Approximately  half  the  solar  irradiance  occurs  in  the  infrared  portion  of  the  spectrum  and  most  of  this 
is  absorbed  and  converted  to  heat  in  the  upper  few  centimeters.  The  remaining  half  of  the  incident  solar 
irradiance  occurs  in  the  visible  and  ultra-violet  port.on  of  the  spectrum  The  flux  of  solar  irradiance  that 
penetrates  below  the  ocean  surface  is  regulated  by  optical  absorption  by  the  water,  phytoplankton,  and  suspended 
particles.  This  absorbed  irradiance  is  mostly  converted  to  heat.  Numerical  modek  have  shown  the  importance 
of  direct  subsurface  solar  heating  in  the  upper  ocean  on  vertical  mixing  Vertical  heating  of  the  water  coliunn 
has  a  significant  impact  on  its  stability.  Enhanced  erosior  of  the  seasonal  thermocline  can  occur  in  areas  where 
substantial  subsurface  heating  occurs,  as  has  been  shown  for  the  Pacific  The  depth  dependence  of  solar 
heating  is  dependent  on  the  wavelength  of  the  incident  solar  irradiance  and  the  spectral  diffuse  attenuation 
coefficient  of  the  water 


In  open  ocean  water,  the  diffuse  attenuation  coefficient  is  controlled  by  the  phytoplankton  concentration 
Morel  ^  has  derived  relations  between  chlorophyll  concentration  and  the  diffuse  attenuation  coefficient  in  the 
visible  portion  of  the  solar  spectrum  from  in  situ  measurements.  The  spatial  variability  of  the  phytoplankton  has 
been  shown  to  be  patchy,  and  the  temporal  distribution  of  the  phytoplankton  has  been  shown  be  h^hly  variable 
based  on  satellite  ocean  color  observations  The  solar  irradiance  distribution  is  ako  highly  variable.  Thk 
variability  k  primarily  a  consequence  of  variations  in  cloudiness  and  aerosol  concentrations  ^  The  patchiness  of 
both  solar  irradiance  and  phytoplankton,  causes  spatial  variations  in  solar  heating. 
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The  objective  of  this  effort  is  to  define  the  visible  irradiance  incident  at  the  sea  surface  and  how  it  is  absorbed 
witlun  the  surface  waters  of  the  North  Atlantic.  This  paper  will  define  the  contribution  to  heating  of  the  water 
column  for  various  visible  wavelengths  resulting  from  the  vertical  spectral  distribution  of  irradiance.  The  North 
Atlantic  in  May,  1979  was  selected  for  analysis  because  of  the  elevated  phytoplankton  concentration  associated 
with  the  spring  bloom.  High  solar  irradiance  is  also  observed  at  this  time  of  year  in  the  North  Atlantic. 


The  general  form  of  the  equation  for  the  solar  heating  rate  of  the  upper  ocean  is  expressed  as; 

dT  IQ  W, 

dt  pCp 

where  is  the  downward  diffuse  attenuation  coefficient  and  is  the  downward  solar  irradiance  2t  the  sea 
surface.  Both  of  these  parameters  are  dependent  on  wavelength,  p  is  the  water  density,  Cp  is  the  heat  capacity 
of  water,  and  z  the  is  water  depth. 


The  diffuse  attenuation  coefficient  that  defmes  the  rate  of  decay  of  irradiance  in  the  water  column  is  dependent 
on  the  vertical  distribution  of  chlorophyll  in  the  open  ocean  Our  computations  assume  the  chlorophyll 
concentrations  which  are  observed  from  satellite  ocean  color  are  constant  within  the  upper  35  m.  Morel  ^ 
defined  the  spectral  dependence  of  based  on  chlorophyll  as: 

ka^)  -  XcW  c 

where  chi,  and  epsilon  are  constants.  His  results  indicate  that  in  water  of  low  chlorophyll  concentration  the 
minimum  occurs  at  about  490  nm  and  that,  as  the  chlorophyll  concentration  maeases,  the  minimum 
shifts  to  about  580  nm  (Figure  1). 


The  spectral  visible  solar  irradiance  at  the  sea  surface;  Eq^;  is 
controlled  by  the  extraterrestrial  solar  spectrum,  and  the  spectral 
transmittance  by  the  atmosphere,  which  is  primarily  governed  by 
the  Rayleigh  scattering  in  the  atmosphere  and  the  ozone 
distribution  The  magnitude  of  the  incident  solar  irradiance  also 
depends  on  time  of  year,  time  of  day,  and  latitude. 


By  coupling  the  spectral  irradiance  at  the  sea  surface  (EgJ,  and 
the  corresponding  spectral  diffuse  attenuation  coefficient  of  the 
water  column  downward  irradiance  at  any  depth  can 

be  computed  using  Beers  law.  The  heating  rate  can  then  be 
computed  assuming  that  within  the  visible  domain  the  absorbed 
radiance  at  a  depth  is  essentially  totally  converted  to  heat  Heat 
associated  with  photosynthetic  storage  processes  is  very  small  (<2 
%)  and  has  been  neglected 

2.  METHODS 


chlorophyll  concentration  from  Morel 
(1988) 


The  Coastal  Zone  Color  Scanner  (CZCS)  global  processing  data  set  was  used  to  obtain  a  monthly  composite 
of  the  diffuse  attenuation  coefficient  at  490  nm  for  May  1979  The  monthly  optical  properties  of  the 
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North  Atlantic  were  obtained  from  the  water  leaving  radiance  at  443  and  550  nm  from  CZCS  ’  (Figjire  2a).  The 
mean  for  approximately  20  km  squares  (at  the  equator)  were  computed  from  of  all  cloud  free  CZCS  data 
during  May  1979.  There  are  several  central  Atlantic  areas  where  no  CZCS  data  were  obtained.  The  North 
Atlantic  illustrates  high  attenuation  coefficients  assodated  with  the  sprii^  bloom  and  the  North  Sea.  The 
southern  and  central  areas  of  the  Sargasso  Sea  indicate  very  clear  water,  i.e.,  low  K,,  values.  Notice  the  patchiness 
of  the  optical  properties  especially  in  high  K„  areas  such  as  coastal  and  upwelling  regions  (African  coast).  The 
mean  montUy  data  were  converted  to  Kjjdo,  K,j44o,  KdS40>  K<is9o»  K,i64o.  Kdew  using  the  relationships  of 
Austin  and  Petzold^^.  These  visible  wavelengths  were  used  to  represent  heating  rates  in  the  visible  spectrum. 


The  incident  solar  spectral  irradiance  was  computed  for  May  1979  for  the  North  Atlantic  using  an  extension  of 
Bird’s  model^'^  (figure  2b).  This  model  was  used  to  compute  the  hourly  incident  spectral  irradiance  at  seven 
visible  wavelengths  (defined  above)  using  (1)  the  aerosol  concentration  from  La670  of  CZCS,  (2)  the  ozone 
concentration  from  the  Total  Ozone  Mapping  Spectrometer,  and  (3)  cloud  cover  from  the  Real  Time 
Nephanalysis  (RTNEPH)  Air  Force  3-hourly  observations.  These  data  bases  from  May  1979  were  used  to 
compute  the  monthly  average  solar  irradiance  at  locations  coincident  with  the  monthly  composite  diffiise 
attenuation  coefficient.  The  spectral  dependence  of  the  solar  irradiance  model  is  described  in  Terrie  and  Amone 
^  in  this  issue.  The  importance  of  aerosols,  in  addition  to  clouds,  is  emphasized  for  certain  locations  in  the  North 
Atlantic  (e.g.,  the  African  Sahara  coast).  The  irradiance  computed  at  the  sea  surface  is  in  units  of  W/mVmn  and 
the  values  obtained  compare  with  other  models'^  Fi^re  2b  illustrates  elevated  irradiance  at  490  nm  at  20°  N, 
just  west  of  the  African  coast.  A  general  decrease  is  observed  to  the  north.  A  band  of  decreasing  solar  irradiance 
is  observed  along  the  Intertropical  Convergence  Zone  because  of  high  cloud  cover.  Similarly,  notice  the  variability 
in  the  irradiance  distribution  at  scales  less  than  100  km. 


3.  RESULTS 

The  average  monthly  irradiance  was  computed  at  0, 5, 10, 15, 20,  and  35  m  depth  for  seven  wavelengths  (defined 
above).  Figure  3  shows  the  monthly  mean  irradiance  field  at  490  nm  at  depths  of  5  and  35  m.  At  5  m  the  490 
nm  wavelength  has  a  similar  distribution  to  the  surface  irradiance,  especially  in  the  clear  waters  of  the  Sargasso 
Sea.  In  the  more  turbid  areas,  such  as  the  North  Sea,  the  irradiance  field  is  strongly  absorbed  and  the 
distribution  appears  similar  to  that  of  the  distribution.  Large  patches  of  phytoplankton  located  south  of 
Greenland  have  a  significant  affect  on  the  irradiance  field  at  5  m.  Very  high  irradiance  levels  are  observed  as 
a  "hot”  area  at  20°  N  in  the  central  North  Atlantic.  This  area  has  high  surface  irradiance  and  low  water 
attenuation.  Similar  areas  of  high  incident  irradiance  and  low  attenuation  water  have  been  observed  in  the 
Pacific.  In  the  northern  latitudes,  high  values  are  observed  with  decreased  solar  irradiance.  The  subsurface 
irradiance  is  rapidly  attenuated  in  the  surface  waters  and  closely  resembles  the  water  attenuation  distribution. 


The  35-m  irradiance  field  at  490  nm  is  very  different  from  the  surface  distribution.  The  northern  latitudes 
indicate  a  substantial  reduction  in  the  irradiance  distribution  as  a  consequence  of  reduced  irradiance  and 
absorption  by  elevated  chlorophyll  concentrations.  The  turbid  North  Sea  waters  have  little  irradiance  penetrating 
to  35  m.  The  clear  waters  at  20°  N  again  show  the  "hot"  area  suggesting  oligotrophic  conditions.  Near  the  US 
east  coast  the  35'm  irradiance  field  shows  the  mesoscale  features  associated  with  the  Gulf  Stream.  The  35-m  field 
distribution  closely  resembles  the  patterns  of  the  data. 


Comparisons  of  spectral  irradiance  at  10-m  depth  at  490  nm  with  440, 540  and  590  nm  are  shown  in  scattergrams 
for  the  entire  North  Atlantic  (figure  4).  The  histograms  for  the  entire  North  Atlantic  are  also  shown  for  each 
scattergram.  The  490,  540,  and  590  histograms  show  a  unimodel  distribution  centered  at  03,  0.25,  and  0.05 
W/m^/nm  from  the  "hot"  area  at  20°  N.  The  histogram  of  the  440  nm  distribution  is  bi-model.  Note  that  the 
higher  irradiance  distributions  are  associated  with  clear  waters  and/or  low  latitudes,  and  that  lower  irradiance’s 
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represent  coastal,  high  values  or  high  latitude  waters.  The  440  -  490  comparison  shows  that  clear,  high 
irradiance  areas  have  a  nearly  covarying  (1:1)  irradiance  relationship,  whereas  at  low  irradiance  levels  (high  K,j 
waters)  the  440  irradiance  is  strongly  absorbed  and  is  lower  than  490  (05:1).  The  540  to  490  comparison  shows 
a  reverse  relationship  with  low  irradiance  (higher  K,,  waters)  having  slightly  higher  540  irradiance  than  490 
(approximately  a  1.2:1).  Conversely,  similarly  for  high  irradiance  (low  waters),  the  540  irradiance  is  lower  than 
the  490.  In  turbid  water  increased  transmittance  of  light  shi^  to  the  higher  wavelengths;  540  nm  and  is 
responsible  for  the  increase  in  540  irradiance  over  490.  At  590  nm  the  irradiance  is  confined  to  the  surface 
layers.  The  490  -  590  irradiance  comparison  for  the  entire  North  Atlantic  indicates  very  low  irradiance  levels  at 
10  m  and  a  monotonic  relationship  with  the  490  irradiance. 


The  average  heating  rate  in  the  upper  10  m  was  computed  from  the  substu-face  irradiance  field.  The  average 
heating  rate  (deg  C/day/nm)  at  4%  and  590  nm  is  shown  in  Figure  5.  The  490  nm  heating  rates  (figure  5)  are 
lower  in  the  tropical  Atlantic  even  though  the  irradiance  levels  are  shown  to  be  higher  (figure  2b).  The  low  K„ 
of  these  waters  results  in  a  substantial  amount  of  solar  irradiance  penetrating  below  the  upper  10  m.  In  northein 
latitudes,  the  heating  rate  at  490  nm  is  higher,  especially  in  the  high  K„  waters  of  the  North  Sea.  Although  the 
incident  irradiance  generally  is  lower  in  the  north,  note  that  the  heating  rates  in  the  upper  10  m  are  substantially 
higher  than  at  the  lower  latitudes.  The  590  nm  heating  rate  (figure  5b)  is  high  throughout  the  North  Atlantic  (see 
histogram  Figure  4)  showing  little  variation.  The  upper  10  m  absorbs  most  of  the  590  nm  irradiance.  The  smaU 
scale  variability  that  is  observed  in  the  figure  is  related  to  incident  surface  irradiance  and  not  to  the  attenuation 
coefficient  of  the  water. 


Spectral  heating  rates  were  computed  for  the  North  Atlantic  for  depth  increments  of  0-5, 5-10, 10-15, 15-20,  and 
20-35  m  (Figure  6).  The  mean  heating  rate  in  the  upper  5  m  is  greatest  at  640  nm  because  of  the  combination 
of  high  incident  solar  irradiance  at  the  sea  surface  and  high  attenuation  in  the  water  (Figure  1)  at  this 
wavelength.  At  wavelengths  above  640  nm  the  heating  rate  in  the  upper  5  m  is  reduced  because  of  lower  solar 
irradiance.  At  wavelengths  below  640,  the  heating  rate  is  reduced  because  of  the  lower  attenuation  and  the 
resulting  increased  penetration  of  irradiance  below  5  m.  At  the  deep  depths,  the  maximum  heating  rate  shifts 
to  440  -  490  nm  because  of  high  solar  irradiance  and  low  attenuation  in  the  water  within  this  spectral  range.  At 
the  low  to  moderate  chlorophyll  concentrations  typical  of  most  of  the  North  Atlantic,  the  minimum  attenuation 
of  solar  irradiance  occurs  at  450  -  500  nm  (Figure  I).  At  wavelengths  below  390  nm,  the  amount  of  available 
incident  solar  irradiance  decreases  significantly  in  addition  to  the  penetration  depth  of  irradiance. 


4.  SUMMARY 

This  effort  represents  a  convolution  of  the  spectral  distribution  of  solar  irradiance  and  the  spectral  attenuation 
of  ocean  waters.  The  spectral  irradiance  at  different  water  depths  for  the  North  Atlantic  during  May  1979  was 
calculated  using  satellite  derived  solar  irradiance  and  water  diffuse  attenuation  coefficient  (K^).  These  monthly 
climatologies  represent  synoptic  characterization  of  the  North  Atlantic  during  the  spring  bloom  when  chlorophyll 
concentration  are  high  and  the  water  diffuse  attenuation  coefficient  is  elevated.  The  incident  solar  spectral 
irradiance  and  the  water’s  diffuse  attenuation  coefficent  were  used  to  conpute  teh  irradiance  at  5  and  35  m.  The 
irradance  distributions  at  5  m  was  similar  to  the  incident  solar  irradiance  especially  in  the  trypical  Atlantic  where 
the  water  was  small.  A  elevated  irradiance  "hot"  area  was  observed  at  20°  N  where  clear  water  permited 
significant  irradiance  penetration  to  depth  (35  m).  North  latitudes  had  reduced  inciden  irradiance  and  elevated 
attenuation  from  the  high  chlorophyll  blooms.  Minimal  irradiance  was  observed  at  depth  at  these  latitudes. 
Spectral  comparisons  of  the  490  nm  irradiance  and  440  and  590  nm  distribution  at  10  m  were  illustrated  in 
scattergrams  and  histograms.  The  440  and  490  nm  irradiance  fields  are  similar  in  low  K,,  waters  which  comprise 
the  majority  of  the  North  Atlantic.  Differences  occur  in  the  higher  K„  waters  (coastal)  where  the  irradiance 
distribution  are  less  penetrating.  With  the  assumption  that  the  subsurface  irradiance  is  totally  absorbed  and 
converted  to  heat,  the  spectral  heating  rates  were  computed  for  five  depth  increments  (0-5,  5-10,  10-15,  15-20, 
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and  20-35  m)  at  seven  wavelengths  (390, 440, 490, 540,  590, 640,  and  690  nm).  The  spatial  variability  of  the  bio- 
optical  properties  coupled  with  the  spatial  variability  of  the  surface  solar  irradiance  distribution,  results  in  a  wide 
variability  in  the  heating  rates  in  the  northern  and  tropical  North  Atlantic. 


The  general  distributions  observed  in  this  analysis  show  that  major  latitudinal  differences  occur  in  heating  rates 
in  addition  to  small-scale  variability.  The  low  K,,,  high  irradiance  waters  observed  at  20°  N  show  that  substantial 
solar  heating  occurs  at  depth  in  these  regions.  Contrastingly,  the  northerly  waters  with  high  K,,,  show  that  very 
little  heating  occurs  at  depth. 


Additional  efforts  are  required  to  understand  the  temporal  variability  of  spectral  heating  rates.  Improved  satellite 
sensors  will  provide  new  looks  at  these  problems  to  help  understand  the  coupling  of  air-sea  interactions  and  their 
affects  on  subsurface  thermodynamics. 
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Figure  4.  Comparison  of  the  spectral  irradiance  at  10  meter  depth  for  the  entire  North  Atlantic  is  shown  in  these 
scattergrams.  Coincident  440,  540  and  590  nm  irradiance  values  (y-axis)  are  compared  against  the  490  nm 
irradiance  values  (x-axis)  for  the  entire  image  field.  The  histogram  of  the  North  Atlantic  irradiance  is  shown  in 
on  the  side  of  the  axis.  The  relationship  shown  by  these  curves  is  related  to  the  spectral  attenuation  of  water 
column.  The  scatter  in  the  relationship  is  related  to  differences  in  the  solar  spectral  irradiance  distribution.  The 
peak  of  the  irradiance  observed  in  the  histograms  is  associated  with  ven,  clear  waters  in  the  central  south  Atlantic 
cast  of  Africa  where  a  "hot"  spot  occurs. 
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Spectral  Heating  in  North  Atlantic 


Mean  heating  Rate  within  Depth  Increment 


Figure  6:  The  Mean  Spectral  Heating  Rate  in  the  North  Atlantic 
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